Instantaneous coherent destruction of tunneling and fast quantum state preparation 

for strongly pulsed spin qubits in diamond 
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Qubits driven by resonant strong pulses are studied and a parameter regime is explored in which 
the dynamics can be solved in closed form. Instantaneous coherent destruction of tunneling can be 
seen for longer pulses, whereas shorter pulses allow a fast preparation of the qubit state. Results 
are compared with recent experiments of pulsed nitrogen-vacancy center spin qubits in diamond. 



I. INTRODUCTION 

Quantum systems with well-isolated two-level subsys- 
tems that can be addressed independently, with states 
that can be initialized and manipulated, enjoy much in- 
terest nowadays because of possible use as quantum bits 
in quantum information processing. Quantum-state ma- 
nipulation is usually described by time-dependent Hamil- 
tonians, and the time-dependent (driving) part of the 
Hamiltonian can either be weak or strong. Weak driv- 
ing is the default situation, but here the focus will be on 
strongly driven quantum systems. 

In the weak-driving regime, the effect of faster and 
larger-amplitude driving may be the same as of slow 
small-amplitude driving, as only the area under the pulse 
is important. This is well known for example for 7r-pulses, 
which are used to invert the population of a qubit. Since 
decoherence sets the time scale within which quantum 
information processing should occur, it is then a good 
idea to speed up quantum state manipulation by mak- 
ing pulses shorter, with an amplitude of the pulses that 
grows concomitantly to keep the same pulse area. This 
strategy works until the peak amplitude grows so large 
that the driving can no longer be considered weak, and 
the dynamics becomes more complex. 

Nevertheless, it may be a good idea to go to the strong- 
driving regime to seek new ways to manipulate qubits 
fast compared to the decoherence time [l|. The strong- 
driving limit has not been studied extensively in a quan- 
tum information perspective. Important was the discov- 
ery of dynamical decoupling of open quantum systems 
from their environments [J Q • The recent advent of 
novel types of qubits with strong couplings to cavities and 
external fields naturally drives the interest in strongly 
driven quantum systems, for example in the context of 
amplitude spectroscopy of superconducting qubits [3-0] 
and optomechanical systems [8|. 

Actually, the interest in strongly driven quantum sys- 
tems dates back before quantum information theory. The 
strong-coupling regime was studied for example in chemi- 
cal physics, in order to control a molecular reaction by in- 
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tense laser pulses 0, [l(| . Many strongly driven quantum 
systems exhibit common phenomena for which there is 
no analogy in the weak-driving limit. One such universal 
phenomenon is coherent destruction of tunneling (CDT), 
which was discovered theoretically and explained in 1991 
by Grossmann, Dittrich, Jung, and Hanggi [ll|,[l2j]. For a 
qubit, it is the phenomenon that a tunneling between the 
two quantum states is brought to a standstill by strong 
driving, but only for specific values of the quotient of 
driving amplitude and frequency (details given below). 
For other parameter values, the tunneling amplitude is 
renormalized, though not 'destructed', due to the driv- 
ing. An explanation of CDT in terms of destructive inter- 
ference of multiple Landau-Zener-Stiickelberg transitions 
was given by Kayanuma soon after [l3[ , see also [3, [l5j] . 

CDT is currently actively studied in several sub-fields 
of physics. Here I only mention a few recent devel- 
opments. The recent single-particle tunneling exper- 
iment on strongly driven metastable argon in optical 
double-well potentials fl6j is a very direct observation 
of CDT, and close in spirit to the original theoretical pa- 
pers [ll|, [l2| • Very interesting are recent many-bo dy gen- 
eralizations of CDT, both theoretical predictions 17H21 1 
and experimental demonstrations [22, 23]. For example, 
for a driven Bose-Einstein condensate in a double well, 
many-body CDT phenomena are predicted that sensi- 
tively depend on the number of particles that have al- 
ready tunneled [20(. Furthermore, for Bose-Einstein con- 
densates in shaken optical lattices, it was predicted [I?} 
and shown experimentally (22[ that the shaking-induced 
renormalization of the tunneling parameter enables the 
switching between the superfluid and the insulator state. 
Interesting is also the recent proposal for routing of quan- 
tum information by employin g co herent destruction of 
tunneling in chains of qubits [2l[ . Moreover, coherent 
destruction of tunneling is also seen in opt i cs, in particu- 
lar in coupled- waveguide structures [TH, [24- 26] , for which 
the coupled spatial wave equations have the same math- 
ematical form as the temporal equations of motion for 
quantum systems exhibiting CDT 1271. F urther work on 
CDT can be found in the reviews [28l - l30l ]. 

The standard system for observing CDT is a harmon- 
ically driven two-level system. With fast manipulation 
of qubits in mind, here we consider driving by pulses in- 
stead, allowing for several types of pulse envelopes, and 
discuss the possibility of observing and exploiting CDT 
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for quantum information processing in this more com- 
plex case. Related work on pulsed systems can be found 
in 0,0, EE). 

We then apply our analysis to nitrogen- vacancy (NV) 
center spin qubits in diamond strongly driven by mi- 
crowave pulses. This is a very promising type of solid- 
state qubit for several reasons [32] . For example, its 
coherence time is long even at room temperature, and 
its state can be initialized by optical pumping. And 
one can transfer its state to even longer-lived nuclear 
spins [H, Hil . NV center spin qubits may become build- 
ing blocks of a quantum repeater [35[. Achieving coher- 
ent coupling between NV centers and superconducting 
flux qubits is another exciting possibility [36| . 

Although our results are general, the very recent mea- 
surements by Fuchs et al. [l| have shown a great level 
of control of the fast and strong microwave pulses by 
which NV center spin qubits can be prepared in desired 
quantum states, and were the main motivation for car- 
rying out the research presented here. We discuss under 
what circumstances CDT in NV center spin qubits could 
be observed, and present analytical results that can be of 
help to prepare hnal quantum states using short resonant 
pulses. 

The structure of this article is as follows. We introduce 
our model in Sec. [TT] with dynamics and some approx- 
imations. This is followed by the exploration of CDT 
phenomena in strongly pulsed qubits in Sec. lIIII The ap- 
proximate treatments and numerically exact results are 
compared in Sec. IIVI The choice of optimized pulses is 
briefly discussed in Sec. |V] The application of the analy- 
sis to nitrogen- vacancy centers in diamond can be found 
in Sec. I VII before the conclusions in Sec. IVIII 



II. QUBIT DRIVEN BY A STRONG PULSE: 
MODEL, DYNAMICS, AND APPROXIMATIONS 

A. Model and equations of motion 

Consider a two-level system with states {|0),|1)} 
driven by a pulse, as described by the Hamiltonian 



luq = 0, the excited-state population Pi(t) varies between 
and 1, but if A/ujq -C 1, then the excited-state popula- 
tion stays negligible before the driving is switched on. For 
A(t) = A constant, the driving is harmonic, and Eq. (fTJ) 
in the particular case of constant-amplitude driving will 
be referred to as the CDT Hamiltonian in standard form, 
or standard CDT Hamiltonian. With the goal in mind 
of fast interactions with qubits, in the following we con- 
sider instead the qubit dynamics under strong driving 
with short pulses. 

With the help of the time evolution operator Uo(t) = 



exp 



■(i/H)a z fdTh z (T) 



an interaction picture can be 

defined in which the Hamiltonian H{t) is given by 
UQ(t)Aa x Uo(t). In matrix representation for the state 
\ip(t)} = ci(£)|l) + co(£)|0), the coupled equations of mo- 
tion for the two coefficients become 



ih 



A(t) 
A*(t) 



where the dots denote time 
time-dependent interaction is 

A exp 



derivatives and 
defined as A(t) 



(2) 



the 



2i/H f AtH z {t) 



These coupled equations of mo- 
tion can be solved numerically without further ado, but 
some striking features of the dynamics can be understood 
after only a little more analysis. 



B. Rotating-wave approximation 

Driving with a pulse adds some complexity to the dy- 
namics as compared to harmonic driving. To keep the 
presentation as simple as possible, let us first choose a 
very convenient pulse shape, obtain some analytical re- 
sults, and then argue that the results can be generalized 
to other pulses. We choose the simple envelope function 
A(t) — Ap(t), with maximum amplitude A and dimen- 
sionless pulse shape function pit) = exp(— |i|/T p ), a two- 
sided exponential with decay rate r" 1 . In this particular 
case, the integral in the exponents of Eq. ^ exactly be- 
comes 



H(t) = h z (t)a z + h x a x = -[uj + V(t)]cr z + hAcr x 



wo/2 + V(t)/2 
A 



A 

-ujo/2 - V(t)/2 



(1) 



where ujq is the energy difference and A the interac- 
tion strength between the two states of the undriven 
qubit, and V(t) — A(t) cos(wi) is the driving field with 
frequency u> and pulse envelope A(t). The a x>z are 
Pauli operators, with a x defined as |0)(1| + |1)(0| and 
<j z = |1)(1| — |0)(0|. For vanishing driving amplitude or 
as long as A(t) is still identically zero after the initial time 
to, the Hamiltonian (TT]) is constant and describes simple 
tunneling dynamics with sinusoidal population exchange 
between the two levels at a frequency i/^o + 4A 2 . For 



2i/h J dTh z (r) = iuJot+icos(x)A(t)sm(ujt^fx)/ u 'j (3) 

with x = arctan[l/ (wr p )] and =F = — t/\t\. The amplitude 
A is reduced by a factor cosx, and x a l so shows up as a 
phase factor. Our convenient pulse shape enables us to 
proceed as in the analysis of CDT for harmonic driving 
by using the mathematical identity known as the Jacobi- 
Anger expansion, 



tx sin a 



E 



J n (x) e l 



(4) 



where J n denotes as usual a Bessel function of the first 
kind. The upper-right matrix element Hiq in Eq. ^ 
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then becomes 

oo 

A Jn[cos{x)A{t)/uj]e^ mx e l{uJo+n ^ t , (5) 

n— — oo 

where the =p signifies that the phases in exp(=pznx) flip 
sign at t = 0. We have not made approximations yet. We 
will make one now, based on the widely different time de- 
pendencies of the terms in the expansion ([5]). One can 
choose the driving frequency uj such that there is an n- 
photon resonance, i.e. ujq + n Tes uj = for some integer 
rt rcs , and zero- and one-photon resonances (n ros = 0, — 1) 
are most common. For harmonic driving this makes 
the n res -term stationary, while all others oscillate. All 
other terms can be neglected when even the slowest other 
terms, the (n res ±l)-terms, average out on the interaction 
time scale A^ 1 . This is a rotating- wave approximation 
(or RWA) and for harmonic driving it is valid if uj 3> A. 
Throughout this paper we will indeed make this assump- 
tion of high-frequency driving. 

In our case of pulses, the situation is slightly more com- 
plex, since the Bessel-function coefficients have become 
time-dependent, including the resonant n rcs term. For 
slowly-varying pulses, the almost-stationary term will 
vary much slower than all other terms, and it is intu- 
itive that it is still a good approximation to keep only 
this term. For very short pulses containing only a few 
oscillations, the validity of the approximation of keeping 
only one time-dependent term is not so clear, and we 
will test it numerically below. (A thorough discussion 
of the RWA of pulsed systems could be given using an 
adiabatic theorem for Floquet-Bloch states [H,|3(J-) The 
RWA dynamics is described by the equations of motion 

AJ„ ros [co S ( X )A(t)M a x ( £ ) , (6) 

where we got rid of phase factors by defining c\ — 
exp[±m ros x/2]ci and en = exp[=Fm res x/2]co, and then 
leaving out the tildes [37j. As in the standard analysis 
of CDT, we see that due to the strong resonant driving, 
the interaction A is replaced by an effective interaction 
strength A e g = AJ„ ros [cos(x) A(t)/u!]. Since the Bessel 
functions are bounded by unity, the effective interaction 
is always smaller. 

For harmonic driving, coherent destruction of tunnel- 
ing is a vanishing effective interaction, A c g = 0. It occurs 
for those specific values of A/u) for which J nreB (A/u) van- 
ishes, and recall that n Tes was fixed by our choice of res- 
onant driving frequency u. Actually the name coherent 
destruction of tunneling is usually reserved for the case 
^rcs = and wo = 0, where it is the driving u 3> A that 
stops the tunneling that otherwise exists between the two 
undriven degenerate energy levels. Here the name CDT 
will be more generally used for situations where A e g = 0, 
even if for resonant driving with n rcs ^ 0, tunneling is 
also already suppressed without driving since uq A. 

The essential difference for our pulsed driving is 
of course that the effective interaction A c ff(t) = 



AJ„ ros [cos(x) A(t)/u] has become a time-dependent 
quantity. This means that coherent destruction of tun- 
neling will not occur at all times. Rather, as the pulse 
amplitude is varied from zero to a large maximal value 
A u and back, then A c g(t) vanishes at those in- 
stances in time at which the Bessel function vanishes. 
This phenomenon could be called instantaneous coherent 
destruction of tunneling, or ICDT, and will be studied in 
Sec. IIIII Central result of this section is the equations 
of motion ©, obtained for the double-sided exponential 
pulse shape, with the RWA as the only approximation 
needed to get there. 

C. Slowly- varying envelope approximation 

Now let us see whether the results of Sec. Ill Tit which 
were obtained for a particular pulse shape, carry over 
to more general pulses. Recall that we did the integral 
J* dr h z (r) and the resulting Eq. ([7]) that we obtained for 
the two-sided exponentially decaying pulse depended on 
a parameter \ — arctan[l/(wT p )]. Now if there are many 
oscillations per pulse, then ujt p is much larger than unity 
and x will almost vanish. But if we take the limit \ = 
in the exact solution ([7]), then we find, for ut v — ¥ oo, that 

h z (t) = \f^{u Q+ A(t) C0S {ut)}^ ff+ ^W"*) , 

Notice that there is a simpler way of obtaining the ap- 
proximate result on the right-hand side, namely by tak- 
ing the amplitude A(t) out of the time integral on the 
left-hand side as if it were a constant. This is a slowly- 
varying amplitude approximation (or SVAA). For the 
double-sided exponential pulse we will now derive the 
condition under which the RWA plus SVAA approxi- 
mated dynamics is close to the exact dynamics, for a large 
but not infinitely large number of oscillations (ujt p 3> 1). 
We can make a Taylor expansion to first order in x of 
the left-hand side of Eq. (UJ). This gives the approxi- 
mate result on the right-hand side of Eq. ([JJ plus a term 
±A(t) cos(ut)x/u. We require the latter to be small at 
all times, which gives for the double-sided exponential 
pulses the following combined condition for the SVAA 
dynamics to be accurate: 

uJTp 3> 1 and A/uj <C ujt p . (8) 

Clearly, the second condition on the scaled amplitude 
A/uj is only very weak once the first condition is satisfied. 

We will now generalize our results to other pulse 
shapes, by assuming that Eq. ([8]) gives the sufficient con- 
ditions under which the SVAA gives accurate results for 
all smooth pulse shapes A(t), where r p denotes the typi- 
cal duration of the pulse. In Sec. Ill B I above Eq. © it was 
argued qualitatively that the RWA would be valid if the 
pulses are slow enough. For the validity of the SVAA we 
now find the more quantitative conditions ([8]). It will be 
checked in Sec. IIVI below that this condition also makes 
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the RWA (j6j) valid for pulsed qubits, at least in combina- 
tion with the standard RWA condition A/oj -C 1. After 
making both approximations, and by following the same 
reasoning as in Sec. Ill Bl the equations of motion assume 
the strikingly simple form 

ih^}*AJ^[A(t)/w]vJ*y (9) 

valid for uj A, r" 1 , yj A/t p . If the driving field would 
have been taken as proportional to cos{ut + 4>) instead of 
cos(u>t), then this nonzero driving phase <f> at time t = 
would lead to phase factors exp(±m<^>) that can be ab- 
sorbed in a redefinition of co,i in the same way as we dealt 
with the phase factors exp(±inx) in Sec. Ill Bl Therefore, 
the equations of motion © describe the approximated 
dynamics irrespective of the value of the static phase <j>. 
Our simple description would have lost some of its appeal 
if this had not been the case, since static phases are often 
hard to control experimentally. 

The main results of this subsection are the approxi- 
mate equations of motion (0) , in which the scaled ampli- 
tude A/uj should satisfy the weak condition (JSJ, but oth- 
erwise is a free parameter. Thus the regime of strongly 
pulsed resonant driving is ready to be explored. 

D. Analytical solution of the approximated 
dynamics 

The coupled equations of motion © obtained after 
RWA and SVAA allow a solution in closed form: 

cx{t) = -feexp[i$(i)]+fcexp[-t $(*)], (10) 
co(f) = fciexp[i$(f)]+fc 2 exp[-i$(t)], (11) 

where ki t 2 are constants that are fixed by the state of 
the qubit at the initial time to. Furthermore, in the ex- 
ponents appears the dynamical phase factor 

*(f)=A /*drJ nre JA(r)/a;]. (12) 

J t 

Unless stated otherwise, we will assume that the qubit 
starts in its ground state |0), in which case k\ ; 2 = 1/2 
and the excited-state population becomes 

Px(t) = sin 2 [$(«)] for |V(to)> = |0>. (13) 

As we will see, this solution in closed form can be of 
considerable use. For example, in order to create desired 
final qubit states with the help of strong short pulses, 
one can employ Eqs. (|12II13|) to optimize the pulse shapes 
A(t), as will be explored in more depth in Section fVl 

III. INSTANTANEOUS COHERENT 
DESTRUCTION OF TUNNELING OF A 
STRONGLY PULSED QUBIT 

In Sec. Ill Bl it was anticipated that in strongly driven 
pulsed quantum bits a phenomenon called instantaneous 




-0.5 1 ' ■ ' ■ ' — 

-2000 2000 

— s- cut 



FIG. 1: Instantaneous coherent destruction of tunneling 
(ICDT) seen in a resonantly (u = ojo) and strongly driven 
(A/lj — 10) qubit. The interaction is small, A/tj = 0.075, 
and the pulse width (FWHM) t p = 2000a; -1 . Panel (a) shows 
the envelope function of the pulse with linear rise and fall, (b) 
shows the numerically exact population dynamics of the ini- 
tially (at to — — 3000oj _1 ) unpopulated excited level |1), and 
(c) depicts the Bessel function Ji as a function of the scaled 
driving amplitude A(t)/u). A temporarily frozen dynamics of 
the population around the dashed lines in (b) coincides with a 
zero of the Bessel function in (c), and the arrows in (a) point 
to A(t)/uj = 3.83 and 7.02, indeed the first two zeroes of J_i 
(or of Ji). 



coherent destruction of tunneling (ICDT) would occur. 
Here we will explore that further. In Figure [TJ a pulse 
containing many resonant oscillations (lj = u>o) rises lin- 
early from zero to a maximum amplitude in the strong- 
coupling regime {A/uj 3> 1) after which it falls back to 
zero with the same constant but now negative slope. In 
the weak-coupling regime, the oscillations of the qubit 
population would become faster as A{t) grows, but the 
figure shows that this is certainly not the case during 
this stronger pulse. Rather, the rate of the oscillations 
is seen to correlate with the value of J-i[A(t)/ui] (or, 
equivalently, with —Ji[A{t)/u]), which is in accordance 
with Eq. ©. 

In particular, at the times that this Bessel func- 
tion vanishes, the population dynamics is temporarily 
brought to a standstill. This lasts for only short while, 
as the pulse amplitude keeps on changing, and the term 
instantaneous coherent destruction of tunneling indeed 
seems appropriate for the phenomenon seen in Fig. [1] If 
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one zooms in on the dynamics around an ICDT point, 
then one sees small fast oscillations remaining in the nu- 
merically exact dynamics. This suggests that the near- 
stationary term suppresses the effect of the fast oscillat- 
ing terms, unless its amplitude nearly vanishes, in which 
case it becomes no more important than the neglected 
rotating terms. 

For the linearly rising half of the pulse, the 
time integral over h z (t) is easy to do, and a term 
Acos(ut)/[ui(u>T p )] is neglected against A(t) sin(wi)/w 
when making the SVAA. From these formulae we again 
see that this is all right to make the approximations in 
case LUT p ^> 1 and (A/[ui(u)t p )} <C 1, which confirms the 
assumed generality of these conditions that were derived 
in Eq. ((5J) above for the two-sided exponential pulse. 

It is to be expected and numerics confirms that as 
pulses become shorter while keeping the same maximum 
amplitude, instantaneous coherent destruction of tunnel- 
ing may not be so clearly visible anymore as a temporar- 
ily frozen population. However, whether visible in the 
population or not, as long as the approximate descrip- 
tion remains valid, the effective interaction still vanishes 
at those instances when the Bessel function J ntaa van- 
ishes, and instantaneous coherent destruction of tunnel- 
ing occurs. 
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FIG. 2: Population dynamics Pi(t) for a qubit resonantly 
(oj = u>o) driven with a two-sided exponential pulse with 
maximum scaled amplitude A/ui — 10 and FWHM width 
t p = 40w~ 1 , for A/lu = 0.02 (lower curves) and A/lu = 0.05 
(upper curves). The red solid curve depicts the numerically 
exact dynamics of Eq. lf2l. the purple dashed curve shows 
the dynamics after making the RWA of Eq. ((6J, and the blue 
dash-dotted curve represents the dynamics after also making 
the SVAA of Eq. (fTHJ). The latter two curves overlap. The 
dotted curve is the pulse form p(t). 



IV. NUMERICAL COMPARISON OF EXACT 
AND APPROXIMATED DYNAMICS 

A. Validity of RWA for approximation for 
two-sided exponential pulses 

We will now compare the exact dynamics of Eq. ([2]) 
with the approximated dynamics. In Figure [2] a compar- 
ison is made for a qubit driven by two-sided exponential 
pulses, for which we did a two-step approximation: first 
we made the RWA in Eq. ©, and then also the SVAA in 
Eq. The driving is quite strong but the approximated 
resonant dynamics closely follows the exact curves. The 
approximation is least accurate after the pulse, when the 
exact dynamics of the undriven qubit shows oscillations 
where the approximated population does not vary, and as 
expected these oscillations are larger for larger A. As an- 
ticipated in Sec.HH making the SVAA hardly gives rise to 
extra inaccuracy once the RWA has been made, as exem- 
plified by the overlap of the corresponding curves. Notice 
also in the figure that the discontinuity in the derivative 
of the pulse shape function at time t = does not lead to 
inaccuracies of the approximated dynamics after t = 0. 
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FIG. 3: Population dynamics Pi(t) for a qubit resonantly 
(w = uo) driven with a gaussian pulse with maximum scaled 
amplitude A/lu = 4 and FWHM width r p = 40a; _1 , for 
A/lu = 0.02 (lower curves) and A/lu = 0.05 (upper curves). 
The red solid curves depict the numerically exact dynamics of 
Eq. p]l. while the blue dashed curve shows the approximated 
dynamics after making both the RWA and SVAA as in Eq. ([9]) 
and its solution Eq. (| 13f) . The dotted gaussian is the pulse 
form pit). 



B. Validity of RWA and SVAA for gaussian pulses 

One can also compare the exact and the approximate 
dynamics for gaussian pulses and this is exemplified in 
Figure [3] The exact and approximated curves for A/to = 
0.05 are not quite the same anymore, which shows that 



the condition A/oj <C 1 for the RWA to hold must be 
observed rather stringently. For the smaller interaction 
A/uj = 0.02, the agreement between the two curves is 
good, and we stick to that small value for the interaction 
when discussing pulse shaping below. 
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V. PULSE SHAPING FOR QUANTUM STATE 
PREPARATION WITH STRONG AND SHORT 
PULSES 

A. Weak-driving limit and area law 

The approximated final-state population depends on 
<&(£), the integral given in Eq. (fT2|) . This integral over 
a Bessel function with argument A{t)/u> can be simpli- 
fied if the driving is weak, A/lu -C f, because in gen- 
eral J n {x) ~ x n /(2 n n\) for x <C 1, and also J_„(x) = 
(— l) n J n (x). For weak resonant driving with uj = luq, we 
have rites = — lj an d the time integral then becomes an in- 
tegral over the pulse area J t drA(r). In other words, for 
this specific resonance we find that the excitation prob- 
ability at some time is given by the time integral over 
the pulse up to that time. Therefore, an intended final 
state can be engineered by the right choice of the total 
area under the driving pulse, where the precise form of 
the smooth pulse does not matter. For other resonances 
| fixes | 7^ 1) the corresponding Bessel function has a non- 
linear small- argument behavior, and the corresponding 
'laws' for weak driving are not pulse area laws. 

B. Final states obtained by strong driving with a 
gaussian pulse 

Often the goal of driving a quantum system is to en- 
gineer an intended final state P, HI, H|[ , the pure excited 
state for example. If we consider a qubit driven by gaus- 
sian pulses and take the interaction A small to keep 
the approximate treatment accurate - and let us take 
A/lu = 0.02 - then we have two parameters left to vary 
the pulse, namely the maximum amplitude A and the 
pulse width t p . Figure 0] depicts the probability to end 
up in the excited state as a function of the driving am- 
plitude A, for several pulse widths. From the figure it is 
clear that for the smallest pulse width lut v = 20 the prob- 
ability Pi(oo) is small, whatever the driving amplitude. 
In order to reach complete inversion Pi(oo) = I, which 
according to Eq. ([To]) corresponds to $(oo) = n/2 + mn 
with integer m, one needs lut p = 100 or a bit less. Fig. @] 
shows that small A and t p can not be compensated by 
huge driving amplitudes to obtain Pi(oo) = I. For the 
longer pulses, the final probability becomes oscillatory 
as a function of the scaled driving amplitude A/w, and 
usually one would then choose the smallest value of A to 
obtain the fully inverted final state. For shorter pulses, 
larger values of A are needed. For large-amplitude driv- 
ing the relation between A and r p to obtain Pi(oo) = I 
is not precisely that of an area law. One can say that for 
strong driving the integral (fT2j) over the Bessel function 
replaces the area law Q . 

In summary, a careful choice of the area under a 
pulse can be used to flip the state in the weak-driving 
regime for n rcs = — I . Choosing to satisfy the condition 
sin 2 [$(oo)] = 1 with $ given in Eq. (|T2"j) is a more general 




' • i i i " 1 1 1 — •, • 

2 4 6 8 10 

-> A/lu 



FIG. 4: Final excited-state probability Pi(oo) for a qubit res- 
onantly driven (u = ujo) by a gaussian pulse, as a function of 
the scaled maximal pulse amplitude A/u), computed with the 
approximate solution Eq. (|13|l . Shown are curves correspond- 
ing to different pulse widths r p (FWHM): on the right, from 
bottom to top, the curves correspond to lot v — 20, 50, 100, 
and 150. The interaction A equals 0.02ct> in all cases. 

method that allows the design of stronger pulses with ar- 
bitrary resonance order n res that achieve a state flip of 
the qubit in a shorter time. 

VI. APPLICATION: SPIN QUBIT IN DIAMOND 
DRIVEN BY A STRONG MICROWAVE PULSE 

As an important application, here the Hamiltonian for 
a strongly driven NV center spin qubit in diamond will be 
briefly introduced, and two realizations of the standard 
CDT Hamiltonian with NV centers will be discussed, be- 
fore discussing the theoretical results of previous sections 
in the light of the very recent experiments by Fuchs et 
al. 0. 

The Hamiltonian Hq of the undriven S = 1 ground 
state of the NV center in a static magnetic field B is 
given by 

H Q = ASt+g et i B B-S, (14) 

which is the sum of a zero-field splitting and a Zeeman 
term. The anisotropy term AS Z describes a zero-field en- 
ergy splitting of A — 2tt x 2.88 GHz between the m = 
and the m = ±1 states. The anisotropy implies a fixed 
quantization axis along the axis between the nitrogen and 
the vacancy, which will be called the z axis. In the Zee- 
man term, g e is the electron g- factor and \xb the Bohr 
magneton. Magnitude and direction of the static mag- 
netic field B = (B x , By, B z ) are not specified yet. The 
NV center is driven by a strong microwave pulse 

V(t) = V2A(t) cos(wt) n d ■ S. (15) 

The unit vector nd describes which spin components of 
S = S x x + S y y + S z z the driving field couples to and 
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how strongly. The factor y/2 is taken out for later con- 
venience. Without loss of generality, we can assume that 
rid = (sin 9, 0, cos 9), with 9 the angle with the z-axis. 
The sum of Eqs. (fH|) and (|T5|) is the total Hamiltonian 
H(t) = Ho + V(t) for the driven three-level atom. 

Often the NV centers are operated such that one of 
the three levels can be neglected in the dynamics, and 
the remaining two levels constitute a long-lived qubit. 
For example, one can choose a static magnetic field 
B = (0,0, B z ) that lifts the m = ±1 degeneracy by an 
amount AE and a driving field with a frequency w that 
is resonant with the m = <-> to = — 1 transition and an 
amplitude A -C AE. Then, if the system starts in the 
to = level, then typically the population of the m = 1 
level is negligible at all times. This is the situation in the 
recent experiments by Fuchs et al. Notice that the 
energy E m= o — Eq = for all field strengths B z . After 
redefinition of the zero of energy, the effective two-level 
Hamiltonian is 

H(t) =eer z + A(t) cos(wi) [sin 9 cr x +cos0er z ], (16) 

with e = | A — g e [iBB z \/2. In order to bring this Hamil- 
tonian in the standard form of the CDT Hamiltonian, 
we rotate to a new basis where the Pauli matrices are 
(t x ,t v ,t z ) and where the driving pulse couples to r 2 . 
In the new basis, the effective two-level Hamiltonian be- 
comes (Case 1): 

tt(,s _ ( £cos9 + A{t) cos(wt) esa\9 \ 

' ~ \ esin(9 -scos9- A(t)cos(ujt) J 

This Hamiltonian for a driven spin qubit in diamond is 
of the general form of Eq. ([T]), with the transition fre- 
quency uj given by |e| cos9 = |A — g e ^BB z \ cos(#)/2 and 
the interaction A = e sin 9. The factor \/2 in the driving 
amplitude canceled against the normalization factor of 
Sx ■ It is important to notice that the Hamiltonian (fl7|) 
can easily be changed upon variation of external param- 
eters. For example, the relative quantity A/ujq = tan# 
is fully controlled by the angle of the driving interaction 
with the z-axis, and can assume all values. The magni- 
tudes of Wo and A are both changed upon variation of 
the external magnetic field B z , while the strength A of 
the driving field V(t) is independently controlled by the 
input power. 

Instead of driving the NV center off- axis, one could 
also drive it along the z-axis (9 — 0), and apply a small 
static magnetic field B x in the x-direction. As before, a 
larger static B z field is to be applied to split the to — 
±1 levels. In this case, the Hamiltonian after two-level 
approximation is (Case 2): 

/e + A(t)cos(ujt) g^B \ 

V 9eHsB x -e - A(t)cas(ut) J v ; 

Advantage of the second case may be that A = g e (J-BB x 
and wo = e = (A — g e ^BB z ) /2 can be varied indepen- 
dently, each by their own external parameter B x z . The 



Cases 1 and 2 are two ways to realize the standard CDT 
Hamiltonian for driven spin qubits in NV centers in dia- 
mond. 



A. Comparison with recent and proposals for 
future experiments 

Strongly pulsed driving of NV center spin qubits for 
Case 1 was reported very recently in Ref. [l[, and for 
some strong and fast pulses the qubit changed state very 
rapidly, with possible ultrafast quantum state prepara- 
tion applications. The experiments show that a great 
level of control can be achieved of the large-amplitude 
pulses that drive the long-lived NV center. The param- 
eter regime was not the same as discussed here, mainly 
because the driving was chosen to make an angle cos 9 = 
a/1/3 with the z-axis, which according to Eq. (|T7|) means 
a strong static interaction in the CDT Hamiltonian with 
A/ wo = \/2 ■ The approximate treatment adopted here is 
valid in the rather opposite A/w< 1 limit. As a conse- 
quence, the short pulses that we found in this regime to 
invert the population are not quite as short as found in 
the experiments of Ref. [lj . The advantage on the other 
hand of the i/u < 1 regime is that it can be benefi- 
cial to design the right pulses with the help of analytical 
predictions for their effects on a qubit, such as presented 
here. 

Another difference is that the current maximal exper- 
imental value of A/uj is ~ 1.1. According to Fig. HI this 
enables one to invert qubit populations fast with pulses 
with ujT p — 150 and uj — 50A, while some of the phenom- 
ena discussed above require even higher values of A/uj. 
Thus the driving reported in Ref. |l| is strong, but not 
quite as strong as would be needed to observe coherent 
destruction of tunneling. 

It is therefore an interesting question whether (instan- 
taneous) coherent destruction of tunneling could be ob- 
served in NV center spin qubits. As shown in Fig. [TJ the 
smallest value of A/oj for which CDT would occur for the 
"res = — 1 resonance is 3.83. Achieving this is challeng- 
ing and would require making ujq smaller by increasing 
B z somewhat more, and by working with slightly higher 
pump powers at the lower resonant frequency. In Case 
1 it would also require to drive at a small but nonzero 
angle with respect to the z axis, so that tan6> « 1. In 
Case 2, again smaller resonance frequencies due to larger 
magnetic field B z and larger powers would be needed, in 
combination with a small static B^-field. For fixed A/uj, 
there is freedom to choose A and uj such as to optimize 
the validity of the two-level approximation. Similarly, in 
case of pulsed driving with A/oj somewhat increased and 
A /uj reduced as compared to the recent experiments [l| , 
the ICDT phenomena as discussed in this article are pre- 
dicted to show up in strongly driven NV center qubits in 
diamond. 

CDT by harmonic driving can serve as a useful method 
to lock a qubit with degenerate energy levels in a de- 
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sired state, even in the presence of static couplings that 
otherwise would lead to tunneling. For the NV center, 
the static magnetic field B z can be increased from the 
850 G of Ref. 1] to 1023 G (a slightly higher increase 
than what is needed above), so that e vanishes and the 
rn = 0, — 1 levels become degenerate. In Case I, the tun- 
nel coupling also vanishes with e. Unwanted static tun- 
neling terms A no i SC in the Hamiltonian could nevertheless 
lead to tunneling, but driving with parameters such that 
A c ff = A no i se J (A/uj) — can prevent this. The small- 
est driving amplitude for which this occurs is A = 2.40w. 
A first demonstration of quantum state preservation by 
CDT in NVs seems simpler in Case II, where the effec- 
tive tunnel coupling A e s can be tuned both by the static 
magnetic field B x and by the driving. 

Importantly from a quantum information perspective, 
the qubit state to be preserved can be an arbitrary su- 
perposition, the coherence of which would not get lost 
while driving, at least on time scales that noise in the 
strong driving is negligible. This stabilization applica- 
tion of CDT is the opposite of fast manipulation. Of 
course it is equally important to be able to keep a qubit 
in a certain state as it is to be able to change its state, 
and both can be achieved by appropriately chosen strong 
high-frequency driving. 

VII. CONCLUSIONS 

In summary, the dynamics of qubits was studied that 
are resonantly driven by strong pulses. Starting with 
a simple pulse shape that allowed an analytical treat- 
ment, the parameter regime was explored in which the 
rotating-wave and slowly-varying amplitude approxima- 
tions are valid in case of strong resonant driving. In that 
regime, accurate solutions in closed form were presented 
for the strongly driven dynamics. Coherent destruction 
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